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1. Introduction 

Gauge/gravity dualities have proven to be a remarkable new framework to study a large 
class of strongly coupled gauge theories [1, 2]. However, the gauge theories that are currently 
amenable to such holographic analysis are typically very different from real world QCD. Hence 
constructing a holographic model of QCD remains one of the most important challenges 
for this approach. Currently, the most successful proposal is a construction by Sakai and 
Sugimoto [3, 4] based on a configuration of D8- and D8-branes in a D4-brane background. 
While reliable calculations are limited to large N^ and small Nf/N^, many observables seem 
to show a good approximation to real QCD at low energies. 

A key feature of the Sakai-Sugimoto model is that it exhibits the desired non-Abelian 
chiral symmetry U{N{)l x U{Nf)ji, as well as its spontaneous breaking [3, 4]. Of course, in 
real world QCD, the analogous symmetry is only approximate as it is explicitly broken by the 
quark masses. A shortcoming of the D8/D8/D4 model then is that the quarks are precisely 
massless. While various suggestions have been made to introduce quark masses [5, 6, 7, 8, 9], 
there remain technical difficulties in pursuing these proposals in detail. A recent proposal 
which seems easier to study is based on deforming the model with certain nonlocal operators 



[10, 11, 12]. The underlying microscopic field theory is a five-dimensional gauge theory where 
the chiral quarks are localized on separate four-dimensional defects. Since the fermions of 
different chiralities are separated in the five-dimensional spacetime, no simple local mass term 
can be introduced in the UV field theory. However, this spatial separation can be overcome 
by connecting two quark fields with a Wilson line. Hence a natural suggestion is to introduce 
a nonlocal operator to provide a quark mass deformation [10, 11, 12]: 



[^^(x'^,x^ = 0)]''Pexp 
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L/2 



[V'R(x^x4 = L)]6. (1.1) 



As has been extensively studied for closed Wilson lines [13, 14], such a nonlocal operator 
would be dual to an instantonic string worldsheet which extends between the D8-D8 pair. In 
the following, we examine the properties of these operators in some detail. 

An overview of the paper is as follows: in section 2, we review the construction of the 
Sakai-Sugimoto background. In section 3, we consider the nonlocal mass terms introduced in 
[10, 11]. In particular, we examine the affect of the dilaton coupling to the string worldsheet. 
Even though this coupling only appears at higher order in the a' expansion, we find that in 
the D4-brane background it introduces a interesting modification in the renormalization of 
the Wilson line. In section 4, we explicitly calculate the expectation value of these nonlocal 
operators. In the absence of any local fermion bilinears, this expectation value or condensate 
is an order parameter characterizing the chiral symmetry breaking in this holographic model. 
We close in section 5 with a discussion of our results and by making a few observations 
about future directions. Appendix A provides the details of a calculation of the fiuctuation 
determinant of the worldsheet fields. The latter contributes at the same order as the dilaton 
coupling but does not make any further modifications of the renormalization of the nonlocal 
operators. 

2. Review of Sakai-Sugimoto background 

The Sakai-Sugimoto model [3, 4] is based on the throat limit of intersecting D4- and D8- 
branes, summarized by the array 

0123456789 
D4 : X X X X X 
D8:xxxx xxxxx 
D8:xxxx xxxxx 

The world- volume theory of the N^ D4-branes naturally gives rise to a maximally supersym- 
metric U{N^) gauge theory in five dimensions. Following [15], the x^ direction is compactified 
and antiperiodic boundary conditions are imposed on the fermionic fields around this circle. 
In the far infrared, one might expect that the only relevant degrees of freedom arising from 
this D4 world-volume theory correspond to four-dimensional Yang-Mills with gauge group 



SU{N^). Further the intersection of the D4-branes with iVf D8-branes supports chiral fermions 
in the fundamental representations of the gauge group U{N^) and of the U{Nf) flavour sym- 
metry. These fermions propagate in the 3 + 1 dimensions common to both sets of branes. 
Similarly, the intersection with Nf D8-branes produces an analogous set of four-dimensional 
anti-chiral fermions. Hence, the Sakai-Sugimoto model produces a holographic description of 
QCD in the throat limit of the D4-branes. The dual gravity theory in this framework yields 
reliable results for large N^ and strong 't Hooft coupling. Our current understanding of this 
holographic model is limited to the quenched approximation, i.e., N(/N^ -^ 0, in which the 
D8-branes are probes in the supergravity background. 

2.1 D4-brane background 

Here we review the supergravity background, which we refer to as the the D4 soliton (following 
the nomenclature of [16]). This throat geometry for a stack of Nc D4-branes with antiperiodic 
fermions on the x^ circle is the gravitational dual of a confined phase of the U{N^) gauge theory 
[15], as described above. For comparison purposes, we also consider the supersymmetric D4- 
brane throat with fermions that are periodic on the x^ circle. Both solutions can be expressed 
in the form ^ 

3 3 „ 

ds' = ("ly ( - dt^ + 6,, dx'dx^ + f{u){dx^f) + (^ " (^ + u^ dQ^^ (2.2) 

e'^ = 5.(|)' F^ = 37riV,£^^ . (2.3) 

The four noncompact directions of the gauge theory correspond to t = x^ and x* with i = 
1,2,3, while the coordinate x labels the compact direction. The 56789-directions transverse 
to the D4-branes are described by a radial coordinate u and four angles that parameterize 
a unit four-sphere. The SO(5)-invariant line element on this sphere is dfi^ , and the volume 
form is e^. The function f{u) is given by 



f{u) = 1 - % , (2.4) 



u^ 



but the constant itKK = for the supersymmetric background. 

The D4 soliton appears to have a conical singularity at u = Ukk • Regularity requires that 
the period of the compact direction, x^ ~ x^ + 27r r^ , is given by 

47r i?3/2 

With this choice the x^ circle smoothly shrinks to zero size at u = Ukk- Fermionic fields in 
the bulk must be antiperiodic on this circle, reflecting the antiperiodic boundary condition 



^The normalization for F4, is different from what is prevalent in the literature and has been chosen to be 
consistent with the usual IIA action, I ^ 1^ J d^°x^ (e"^'^(i? + 4(V<;/))^) - ^FI) . 



on fermions in the dual gauge theory. Unhke the sohton background, the supersymmetric 
D4-brane geometry with Ukk = exhibits a naked curvature singularity at u = 0. In that 
case there is no restriction on the periodicity of the x^ direction. Further, the dual gauge 
theory is not confining. 

The supergravity solution described above is completely specified by the string cou- 
pling constant, gg, the RR flux quantum (i.e., the number of D4-branes), N^, and the non- 
extremality constant, Ukk- The remaining parameter is a length scale, R, which is given in 
terms of these quantities and the string length, £s, by 

R'^ = TTgsNJ^. (2.6) 

Various combinations of these parameters have direct interpretations in the dual gauge theory. 
The holographic dictionary gives the five-dimensional gauge coupling as g"^ = {2TT)'^gsis, so 
that the five-dimensional 't Hooft coupling is 

A5 = gi N, . (2.7) 

Since these couplings have dimensions of length there is a power-law running of the dimen- 
sionless effective coupling [17] 

g'^„ = gJN^U, (2.8) 

where the energy scale U is related to the radial coordinate u in the D4-brane throat by 
U = u/2'kII [18]. One finds that the scale of Kaluza-Klein excitations of the compactified 
coordinate x^ gives the characteristic mass for glueballs [19] 



"''■=vrf§-^- p-^' 



1 O ^KK 



Below this scale, the low-energy gauge coupling in four dimensions is ^ym = ffi/^Trr^. 

Supergravity provides a good description of physics in the D4-brane background if two 
conditions are met. First, gravity calculations are reliable if the length scale associated with 
spacetime curvatures is small compared to the fundamental string tension. In the D4 soliton 
solution (2.4) the Ricci scalar has a maximum at u ~ Ukk, where curvatures are of order 
(nKK-R^)~^'^- Hence we require 

^^^^ » 1 • (2.10) 

In terms of gauge theory quantities this can be expressed as the condition 

5ym^^c » 1 (2.11) 

so that the restriction to small curvatures corresponds to a large 't Hooft coupling in the 
effective four-dimensional gauge theory. Second, string loop effects are suppressed as long as 
the local string coupling is small: e''' <C 1. The form of the dilaton (2.3) implies that, for finite 
values of the gauge theory parameters, the inequality e"^ <C 1 can only be satisfied over some 



finite range of the coordinate u. Tlie string coupling eventually becomes 0(1) at a value of u 
given by 

Ucr^t - -i ■ (2.12) 

This critical radius naturally becomes large in the limit (7ym -^ 0. Taken together, equations 
(2.10) and (2.12) indicate that the supergravity analysis in the D4-soliton background is 
reliable in precisely the strong-coupling regime of the 't Hooft limit of the four-dimensional 
gauge theory. 

In the strong coupling regime the QCD scale cannot be decoupled from the compactifi- 
cation scale, e.g., in the confining phase described by the D4 soliton, the QCD string tension 
is T ~ gy^^N^/r^"^ [15, 20]. This means that, for most practical purposes, calculations in 
the holographic framework are reliable in a regime corresponding to a five-dimensional gauge 
theory. Since this theory is nonrenormalizable it should be thought of as being defined with a 
cut-off scale, Uoo = Uoo/2iTa' . Above this scale a UV completion with new degrees of freedom 
is required. This completion may be a lift to M-theory, with the x^^ circle opening up to re- 
veal an asymptotically AdSy x S'^ background (with identifications). On the field theory side 
of the duality, the UV completion of the five-dimensional Yang-Mills theory is given by the 
six-dimensional (2, 0) theory compactified on a circle. An alternative UV completion would 
simply be type IIA superstring theory in the asymptotically flat D4-brane background. 

2.2 D8-brane embeddings 

Our current understanding of the holographic model described in the previous section is 
largely limited to the quenched approximation: N{/Nc -^ 0. In this limit the D8-branes can 
be treated as probes embedded in the supergravity background generated by the D4-branes ^. 
The D8-brane fills the noncompact 0123 directions as well as the angles on the S"^ transverse 
to the D4-branes. The nontrivial aspect of the embedding is given by a function x'^{u) that 
characterizes the D8-brane's profile in the u-x'^ plane. With this choice of embedding, the 
action for the D8-branes is 



I^,^-N,T^,jduu''Jf{u){duX^f+{^ j^. (2.13) 

The resulting equation of motion for x^ (u) is 

^'^^^^^""' ] = . (2.14) 

^f{u){duX^)^ + {R/u)yf{u)J 

The expression within the parentheses is constant. If we assume that the profile is symmetric 
across the u-axis, crosses this axis at some value uq, and is smooth in the vicinity of this point, 
then this constant is given by ^oV^/(^o)- The embedding equation can then be expressed as 

u3/2/(t.)v'<'''/(«)-KS!/(i'0) 



^See [21], for attempts to account for the gravitational back-reaction of the D8-branes. 



The boundary conditions for the D8-brane profile are: asymptotically as u ^ oo, x^ -^ L/2 
and (9mX^ oc l/u^^'"^ -^ 0; at the minimum u -^ uq, x* -^ 0+ and duX^ oc l/(ti — uqY'"^ -^ co. 
The full embedding consists of two halves of this form. Hence the D8- and D8-branes are 
joined in a smooth profile at the minimum radius uq and the two defects are separated by the 
asymptotic distance L in the x^ direction. The limit uq —>■ Ukk yields the trivial embedding 
x^=constant, which in the D4 soliton background corresponds to a smooth joining of the 
D8- and D8-branes with asymptotic separation L = irr^. In the supersymmetric background 
(^KK = 0) the trivial embedding is also allowed with an arbitrary separation L. In this case, 
the D8- and D8-branes terminate on the singularity at u = 0. 

We can gain some intuition for these embeddings by considering the supersymmetric 
background. In this case with u^k = 0, f{u) = 1 and hence the embedding equation reduces 
to 



dx^ i?3/2 



Integrating gives a solution in terms of an incomplete Beta function [22] 



L= . (2.16) 



4, , i?3/2^ r(9/16) i?3/2 /„8 

8^ r(17/16) 8^ 



X [u 



/3(^, 9/16, 1/2) . (2.17) 



This has a finite u ^ 00 limit, so that the asymptotic (coordinate) separation of the D8/D8 

pair is given by 

^3/2 ^3/2 

L = Jim 2x^n) = Co ^j^ ^ 0.7245 -^ , (2.18) 



v/^r(9/16) 



where Cq = 4rd7/i6l • ^^ course, L corresponds to the separation of the defects in the dual 
gauge theory. With the D4 soliton background there is a maximum separation corresponding 
to defects located at antipodes on the x^ circle. In the supersymmetric background the 
x^ coordinate need not be compact, so there is no restriction of this sort. However, other 
considerations bound the maximum value of L which one might consider in this case. Notice 
that increasing L corresponds to smaller values of uq. If uq becomes too small, the D8-brane 
extends into a region of high curvature and the calculation described above is no longer 
reliable. Therefore one can only reliably work with values of L where the minimum of the 
brane embedding is safely outside of this region. 

With periodic fermions, the adjoint sector of the theory is supersymmetric and the gauge 
theory is not confining. In the dual background, free "constituent" quarks are realized as 
strings stretching from uq, the minimal radius of the branes, down to u = 0. The energy of 
these strings corresponds to the mass of the constituent quarks: m^ = UQ/2TTa' . With (2.18), 
this dynamically generated mass scale is given by 

Since 27rr^ > L we always have niq ^ 1/r^ in the strong coupling regime, which reflects the 
fact that the infrared dynamics does not decouple from the compactification scale. In the 



confining background of the D4 soliton background, tliere are no free quarks but one can still 
show that Uq has the interpretation of roughly determining the constituent quark mass niq 
as above [23], at least when no is sufficiently larger than u^k- 

The constituent quarks above are complicated bound states of "current" quarks (i.e., the 
fundamental fields in the UV Lagrangian) and adjoint fields. This is shown in a striking way 
by comparing the quantum numbers of the constituent quarks to those of the current quarks. 
In particular, the current quarks are singlets under the global SO{5) symmetry [3]. However, 
at Uq , the D8-brane is wrapping the internal S"^ and the strings stretching from here to u = 
in the supersymmetric background, can rotate in this internal space. Hence the constituent 
quarks transform nontrivially under 5*0(5). Furthermore, it is likely that quantizing these 
strings will give a spectrum of both bosonic and fermionic states. 

One can also look at the low-lying meson spectrum by considering excitations of the 
world- volume fields on the D8-branes [24, 25]. One finds that this spectrum is (as expected) 
characterized by the mass scale mq/g^g [26], in accord with the standard supergravity formula 
[18]. Here, one explicitly finds that mesons are both fermions and bosons [3], rather than 
just bosons. The latter reflects the fact that these infrared excitations are again complicated 
bound states of both the (fermionic) quarks and (both fermionic and bosonic) adjoint fields 
found in the UV Lagrangian. 

3. Nonlocal mass term 

Recall that the underlying microscopic field theory is a five-dimensional gauge theory with 
(chiral) fundamental matter fields localized on two four-dimensional defects. These defects are 
separated along the x^ circle, so that fermions of different chiralities live at different places in 
the spacetime. Hence a naive mass term of the form [^j^(x^, x^ = — L/2)]"[^r(x'^, x^ = L/2)]^ 
is not possible — in particular, it is not gauge-invariant. As described in the introduction, 
the best one can do is to construct a nonlocal but gauge-invariant operator (1.1) with a 
Wilson line connecting the quarks on the separated defects [10, 11, 12]. This suggests that 
one consider the gravity/string dual as an instantonic Euclidean string worldsheet which sits 
at x^ and extends between the D8-branes in the x^ direction [13, 14, 27, 28, 29].^ Of course, 
this worldsheet does not quite reproduce the operator given above in (1.1). Rather this 
holographic construction introduces an "enhanced" Wilson line which sources both the gauge 
field and the adjoint scalars <1> of the five dimensional gauge theory [13] 



©(x^) = [V'l(x'',x^ = 0)]"Pexp 



l-L/2 

/ {A^-ini^^) dx^ 

J-L/2 



[Vi^(x^x4 = L)]fe. (3.1) 



Here the (constant) normal vector nj indicates the position of the worldsheet in the internal 
space. In principle, one could consider an elaborate contour for the Wilson line connecting 



^A similar class of worldsheet instantons were studied with regard to the U{1)a problem in the Sakai- 
Sugimoto model [30]. 



the two fundamental fermions. However, for the sake of simplicity, we will only consider a 
straight contour (with fixed x^) in the following. 

In analogy with the usual holographic calculations of Wilson loops, the expectation value 
of (3.1) is given by 

(O) ~ e-^ws ^ (3.2) 

where /ws is the worldsheet action for a string stretched between the -D8/D8 pair. However, 
there is an interesting difference between the present calculations and those for the conformal 
super- Yang-Mills theory [13, 14, 27, 28]. As described above, the five-dimensional gauge 
theory under consideration here is defined with a cut-off scale, C/qo = Uoo/2TTa' . The necessity 
of this cut-off is reflected in the nontrivial dilaton profile reaching strong coupling at a large 
radius in the dual D4-brane background. However, the standard Wilson line calculations 
are unaware of this aspect of the physics since those calculations only consider the leading- 
order Polyakov term for the worldsheet action in (3.2). The effect of the nontrivial dilaton 
is first seen at next-to-leading order in a', through the Fradkin-Tsyetlin term [31]. Hence to 
better understand the physics of Wilson line operators (3.1), we are motivated to carry the 
worldsheet calculations to first order in the a' expansion. Consistency at this order demands 
that we also include the fluctuation determinants of the worldsheet fields. 

In the following, we find that the Fradkin-Tsyetlin contribution generates new diver- 
gences in the calculation of (3.2) and hence the renormalization procedure for the worldsheet 
action must be modified. In particular, the new divergences do not seem to be removed by 
the Legendre transform introduced in [27]. We illustrate this in the following sections by eval- 
uating the action for a specific worldsheet in the supersymmetric background. This simple 
calculation exhibits the full set of divergences that appear in subsequent calculations and so 
allows us to give a prescription for renormalizing the worldsheet action to first order in a' . 

3.1 Rectangular Worldsheet 

To get a feeling for the issues that arise in the calculation of (3.2), we consider the trivial D8- 
brane embedding, x^ = ±L/2, in the supersymmetric D4-brane background. This corresponds 
to having the D8- and D8-branes extend straight down along the u direction, from the cut-off 
at Uoo to the curvature singularity at u = 0. The gravity approximation breaks down for 
small values of u, so we introduce an (arbitrary) IR cut-off at u = Uir, with B?''^Uil ^ a' . 
For convenience, we assume that this is implemented by introducing a probe D4-brane where 
the strings can end. The result is the rectangular worldsheet shown in figure 1. 
Let us begin by considering the Polyakov action 

/p = ^ y rf VV5 5"' Gij daX'da' , (3.3) 

where A^ is a worldsheet with boundary dM, gab is the spacetime metric on ^A, and Gij is 
the metric (2.2). Two of the worldsheet scalars can be identified with the coordinates u and 
x^ — the remaining scalars can be ignored for the moment. The worldsheet metric is taken 
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Figure 1: The rectangular worldsheet described in tlie main text. 



to be the same as the pullback of (2.2) to Ai, which is given by 

3/2 



g^.dx^dx" = [^f\dxr + [^J 



du . 



Evaluating the Polyakov action yields 



L 



2'Kl'i 



[Uoo "^IRJ 



(3.4) 



(3.5) 



The first contribution is a UV divergence proportional to the cut-off scale Mqo- This term 
is removed by the Legendre transformation described in [27], which gives a 'renormalized' 
action 
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(3.6) 



Alternately, subtracting the UV-divergent term from (3.5) can be interpreted as a renormal- 
ization of the field theory operator (3.1). 

The Polyakov action is the leading order contribution to the worldsheet action in the a' 
expansion. We must also take into account terms at the next order in this expansion that 
couple to the nontrivial dilaton of the D4-brane background. Specifically, we have to evaluate 
the Fradkin-Tsyetlin term [31, 32] 






(3.7) 



where 7^ is the worldsheet Ricci scalar, s is the proper distance along the boundary dM, and 
/C is the geodesic curvature of the boundary. The latter is defined as 



/C = -eubVat^ 



(3.8) 



where t" and n" are unit vectors tangent and normal to the boundary, respectively. The last 
term in (3.7) is a sum over corners where the embedding of the boundary is not smooth. A 
corner that makes an angle 6 gives a contribution proportional to it — 9, times the value of 
the dilaton at that point. These corner terms can be thought of as arising from 5-function 
contributions to the geodesic curvature. Of course, with a constant dilaton <J>o the sum of all 
the terms in (3.7) gives Ipx = X^O) where x is the Euler character of the worldsheet. The 
worldsheets that we consider all have the topology of a disk, i.e., x = 1- Finally, consistency 
requires that we take into account the fluctuation determinant on the worldsheet at this order 
in a' . This calculation is performed in appendix A, where we find that it does not make a 
significant contribution to the worldsheet action. In particular, these one-loop determinants 
do not generate any additional UV divergences. 

We now evaluate the individual terms in (3.7), beginning with the scalar curvature term. 
The Ricci scalar for the worldsheet metric (3.4) is given by 



ng) 
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4 /?3/2 ^1/2 • 



(3.9) 



and the first term in (3.7) is 
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Next we consider the contributions to (3.7) from the smooth components of the boundary. 
The component of the boundary extending from {—L/2,Uoo) to (L/2,Uoo) has tangent and 
normal vectors given by 



t" 



{R/uoof/\d/dx^r , n'^ = {u^/Rfl\d/du^r 



(3.10) 



Using these expressions in equation (3.8) gives the geodesic curvature along this part of the 
boundary 



/C 



AR \u, 
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1/4 



(3.11) 



The proper distance along this edge is ds = {uoo/R)^'^dx^, so the contribution to the action 
is 
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ds K,^ = - — ^ \—fr] log 



SttR \ R 



9s 



noo\3/4 



R 



(3.12) 
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The component of the boundary between (— L/2, uq) to {L/2, uq) makes a similar contribution; 



it differs by an overall minus sign and the substitution u^ 



The geodesic curvature 



vanishes for the edges of the worldsheet along x^ = it L/2, so they do not contribute to the 
action. Finally we consider the contribution of the four corners of the worldsheet, each of 
which makes an angle 9i = -k 12. Their contribution to the action is 



(3.13) 
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Collecting these terms, the action (3.7) yields 

-'FT rr-r 
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+ log5's 



(3.14) 



Hence the inclusion of (3.7) leads to two new UV-divergent terms in the worldsheet action, 
proportional to ^/u^ and logUoo- These are in addition to the divergent term coming from 
the Polyakov action (3.5). Notice as well that both (3.10) and (3.12) contained potentially 
divergent terms of the form Uoo^'^ log Woo however, these terms cancel out in the final expres- 
sion. 

Although we have used a particularly simple background and worldsheet configuration in 
the present calculation, the structure of the UV divergences depends only on the asymptotic 
behaviour. This means that the result obtained here is in fact universal, and the divergences 
we have found also appear in more general situations. The calculations in section 4 — both 
analytical and numerical — show this explicitly. Therefore, as we discuss below, the results 
of this section lead to a general prescription for renormalizing the worldsheet action and 
obtaining a finite expectation value (O). 

As a final comment here, we note the term log^f^ which arises in (3.14) from the inclusion 
of the Fradkin-Tseytlin term (3.7) in our calculation. Keeping the background scale R fixed, 
(2.6) gives gs oc 1/N^ and hence one finds (O) ~ N^. Given that O is a bilinear of fields in 
the fundamental representation of the gauge group, this latter factor is precisely the expected 
result by the standard large N^ counting. Of course, this factor is a universal result for all 
such worldsheet calculations, as we will see with the examples calculated in the section 4. 

3.2 Renormalization of the Worldsheet Action 

One can try to address the UV divergences in (3.14) by applying the Legendre transform 
described in [27]. The authors there demonstrated that the 'correct' action for observables 
related to the minimal area of a string worldsheet is the Legendre transform of (3.3) with 
respect to some of the loop variables — see also [28]. This is because some of the worldsheet 
scalars satisfy Neumann boundary conditions asymptotically rather than Dirichlet boundary 
conditions. Indeed, as we commented above, implementing this Legendre transformation 
removes the UV-divergent term from the Polyakov action (3.5). However, a straightforward 
application of the same Legendre transform does not cancel the divergent terms in (3.14). 
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To see that this is the case, first vary the full worldsheet action with respect to the 
worldsheet fields. This gives an expression of the form 



51 



d^x^ S^'^Sg^^ + SiSX^ + dxVh ^T^''' 5h^^ + t:i 5X^ 
M L -I JdM L 



(3.15) 



The coefficients of 5g^j_y and 5X^ in the integral over J\A are the worldsheet equations of 
motion, while the coefficients of 5h^y and 5X^ in the boundary integral are the momenta 
pulled back to ^^A. The vr/ are given by 



TT/ 



27r^2 



Gun'^d^X'^ + i^K:di^ 



1 

2^' 



(3.16) 



where n^ is an outward pointing unit vector normal to 9A^, and all fields are evaluated 
at dM. The Legendre transform of the action with respect to some subset {X^^ of the 
worldsheet scalars is denoted / and is given by 



/ = /- f dxVh Y^ TTjX^ . 
JdM ^^jy 



(3.17) 



Following [27], we construct the Legendre transform of Irect with respect to the worldsheet 
scalar u at Uoo 

Irect = Irect " / dx^TTuU . (3.18) 

JdMoa 



Using (3.16), we have 



1 



9 



1 



^" 27r^2 [uj +32 7r i? 3/4^5/4 



The induced metric at Uoo yields v7i = 

dx Vh vr,, u 



[Ur 



,/Rf/^ and so 



L 



dMc 



2vr£2 



Uoo + 



9L 



32 7ri?3/4 



(3.19) 



(3.20) 



The first term cancels the leading power-law divergence coming from the Polyakov term (3.5), 
but the second term does not cancel the corresponding Uqo^ term in (3.14) 



L 



'^rect 



27rj 



■uo + 



9L 



32 7ri?3/2 



/^ - 2 ^) + - ($oo + <I'o) 



(3.21) 



Thus the usual Legendre transform of the worldsheet action does not address the divergent 
terms at next-to- leading order in a'. As the structure of the UV divergences is universal, this 
approach also fails for more general curved embeddings, such as those that we study in the 
next section. 

The simplest method for dealing with the divergences is to subtract the terms in (3.14) 
that depend on Uqo- This approach is closer in spirit to that applied in the holographic 
renormalization of probe D-brane calculations, e.g., see [33, 34]. We have already noted that 
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subtracting the Uoo term in the Polyakov action can be interpreted in the field theory as a 
UV renorniahzation of O and the same interpretation applies to the new terms at next-to- 
leading order in the worldsheet action. Such a subtraction is straightforward for the lioo^ 
term, however, we also have to deal with the logarithmic contribution from the two corners 
at u = Uoo, x^ = ibL/2. An ambiguity naturally arises here because the subtraction, which 
takes the form 3/4 \og{u^ / u^^^) ^ requires the introduction of a subtraction scale u^^^,. Thus, 
our proposal for renormalizing the worldsheet action is 



r(ren) j L 9L Uoo 3 (Uoo\ ,„ ^„x 

The UV divergences that we subtract are universal and render the worldsheet action finite up 
to terms of order a' . As shown in appendix A, there are no divergences associated with the 
fluctuation determinant. In the next section we apply this renorniahzation to the worldsheet 
action for a string stretching between the D8- and D8-branes with the curved embedding. 

4. Worldsheet for smooth D8-D8 pair embedding 

Now we turn to an explicit calculation of the expectation value (O) for the Sakai-Sugimoto 
model [3, 4]. As this holographic model does not permit the construction of a local fermion 
bilinear, this expectation value is the condensate which characterizes the spontaneous break- 
ing of chiral symmetry in this model. As described in section 2, this spontaneous symmetry 
breaking is realized in the gravitational dual by the D8-D8 pair joining together to form a 
smooth U-shaped embedding, as illustrated in figure 2. In the figure, the light blue region is 
the (Euclidean) worldsheet of a string stretching between the D8-D8 pair. The boundary of 
this worldsheet consists of two smooth components: the segment dMoo defined by the cut-off 
Uoo, and the segment dMos defined by the embedding x^ = x^{u). 

Recall that embedding equation (2.15) determines the D8-brane profile x^(n) 

dx^ ^ R^l' u^ TTM 1 .4 ^. 

du u^/'fiu) ^u^fiu)-u^fiuo) • ^' ' 

where uq is the minimum radius where the D8-D8 pair joins smoothly — see figure 2 — and 
f{u) = 1 — (ukk/u) . With Ukk 7^ 0, this equation cannot be solved analytically and so the 
embedding must be determined numerically. 

To numerically solve for the embedding we define the following dimensionless variables 

u ukk , [m x^ 

z = — uj = ^=./ 4.2) 

Uq Uq \ H H 

with Uq the point on the n-axis where the embedding reaches its minimum value. The re- 
striction no < u < Uoo implies 1 < z < Zoo- Similarly, the parameter uj takes values to £ [0, 1]. 
Here the lower limit corresponds to uq ^ 00 but this limit is also realized in the extremal 
background with u^k = 0. The upper bound is reached when the embedding reaches the 
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dM. 




Figure 2: Worldsheet for a smooth embedding of the D8/D8 branes. 



mininiuin radius of the background at u 
the embedding equation becomes 



Ukk- In terms of these dimensionless variables, 



dil){z,uj) z^l"^ Vl - ^^ 



di 



UJ'' 



^ 



3 ^5 



U)-^ z 



l + w^ 



(4.3) 



Solving this equation numerically yields a family of embeddings ip{z,ij) parameterized by lv. 
With the standard dictionary U = u/27t£'^, the parameter to becomes a ratio of scales 
io = -#^ where Ukk and Uq can be thought of as the confinement and chiral symmetry 
breaking scales, respectively. In this model, these are both dynamically generated scales 
determined by the fundamental gauge theory parameters. For example, using (2.5) and the 



subsequent formulae in section 2.1, U^ 



As 



and in the super symmetric background. 



Uq = rUq in (2.19). Hence, in principle lo is also a function of the parameters L, r^ and A5. 
In fact, a relatively simple expression can be derived by first noting that for a sufficiently 
large cut-off Zoo = Uoo/uq, 'ip{zoo,uj) is essentially only a function of u>. Then the asymptotic 
separation L = 2x'^{uoo) can be expressed in terms of '0oo(w) with 



L = 2V'o 




(4.4) 



14 




Figure 3: The function -y/cJ V'oo('j-') that appears on the right-hand-side of equation (4.5). The dashed 
red hne corresponds to L = irr^, where the defects sit at antipodal points on the x^ circle. The dashed 
black line at w = 1 is the bound u = Ukk- 



Using the various expressions in section 2.1, we then find 

-— = Vw^oo(w) 

3 r. 



(4.5) 



and so in fact uj is independent of the coupling A5. Hence the coupling dependence of the 
individual scales C/kk and Uq has canceled in the ratio defining uj. The function y/ujipooi^) is 
shown in figure 3 on the range uj £ [0, 1]. 

Holding L fixed in (4.5) implicitly gives a; as a function of r^^. This allows us to interpret 
ip{z^uj) as a family of D8/D8 embeddings with constant asymptotic separation in Z)4-brane 
backgrounds with different compactification radii r^. Three such embeddings are shown in 
figure 4. Alternately, fixing r^ in (4.5) implicitly gives a; as a function of L. In that case 
^{z,uj) represents a family of D8/D8 embeddings with varying asymptotic separation in a 
fixed, non-extremal Z?4-brane background. 

4.1 Worldsheet Action for the Curved D8-D8 Embedding 

Next we explicitly evaluate the worldsheet action for a string stretched between D8- and 
D8-branes for the curved embeddings described in the previous section. In the extremal 
background the calculation can be performed analytically; in the non-extremal case we must 
use numerical methods. While its contributions are subdominant, we include the Fradkin- 
Tseytlin term in the following for illustrative purposes. 

As in the previous calculation, we identify the coordinates of the Euclidean worldsheet 
with the spacetime coordinates u and x^. The worldsheet metric and dilaton are given by 



ds^-{^f'fiu)id.r^(^y^' 



du' 



uJ f{u) 



^ = log (5s (|)'^ 



(4.6) 
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uo = 179 R 
Ukk = 168 R 



I- Uo = 93 i? 

Ukk = 73 R 

Uo = 53 R 



* x^{u) 



Figure 4: D8/D8 embeddings with fixed asymptotic separation L = R/10, for mkk = (blue), 
Ukk = 73R (violet), and Mkk = 168 i? (green). The non-zero values of Ukk are indicated by red 
dashed lines, and the minimum value of u reached by each embedding is indicated by a black dashed 
line. 



To compute the Fradkin-Tseytlin part of the action we need expressions for the scalar curva- 
ture on A1, as well as the proper distance and geodesic curvature on dAi. The Ricci scalar 
for the metric (4.6) is 

15 ui" 



4^3/2^1/2 ^4^3/2^7/2 " ^ ' ^ 

The normal vector, tangent vector, and geodesic curvature for the component of the boundary 
at u = Uoo are given by 

3/4 



v=(^)"*^J=V V=(-) vTKJV 



/c 



3 (Uoo^ + U 



^ ) 



4i?3/4^_L3/4y7(^ 

On the component of dM. described by the embedding x'^{u) these quantities are 



n,A 






RJ w* 



Ry(uoY^uo) 



u / 



/w 



Moj 



t-u 



i^y^ V^^/N-ng/(no) 



u/ 



nvw 



(4.8) 
(4.9) 

(4.10) 
(4.11) 
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13noV/(no) 
^ = 4i?3/4,i7/4 • (4-12) 

Above in (4.10) and (4.11), the upper (lower) sign corresponds to the portion of the boundary 
with x^ > {x^ < 0). Using these expressions and the diniensionless variables (4.2), the 
worldsheet action is 

i?2 /^l r^- 3 r^ ^ ,, J 3 15a;3\ 

9 , . ^ 1/2 Z^-, ^^\ , 39 Z-^-, yr^7J3 

loTT V z^ y 167r Ji z Vz^ - UJ^Z^ -1+UJ^ 

+^logzoo + log(5.(^)'^') . (4.13) 

This expression must be renormalized according to the prescription in section 3.2, which in 
terms of the diniensionless variables becomes 

(_)_ _j^^ 9L^^ 3 

As described above, this prescription requires choosing a subtraction scale n^ub = uq z^uh- For 
simpUcity, we choose n^ub = uq {i.e., z^ub = 1) in the fohowing. To proceed, we must use (4.4) 
to simphfy various factors, e.g., 



Luo As V4('^) R [uq a 



5 



27r£2 2ttL vr ' H \ R 2tt L 



V'oo(c^) . (4.15) 



In particular, with these expressions, the final result is expressed as a function of the ratio 
A5/L and the parameter lo. Then using (4.14) to explicitly remove the divergent terms from 
(4.13), the renormalized action becomes 

r^-^ - ^ ^P^iu) F,iu) - F,{u) - log (vriV,) + ^ log f %M bf] . (4.I6) 



'WS 



where 



2ttL -^^-^ ' --^ ' ■^^ ' °^ '^' ■ 2 ° V 27r L, 

Fi{uj) = - dz (^00 (^) - i'iz^uj)) + -V'oo(w) (4.17) 

■K Ji \ / vr 

^2(^) = -^ Hdz (^oo(^)-V'(^,^)) ^ + ^Voo(a;) (4.18) 

32 7r Ji V / yjz 8tt 

9uj^ .logzoo 45lj3 r^- logz 

39 Z""^-, Vl-w^ 

dz = log z . 



16 vr Ji z Vz^ - z^uj^ -l+uj^ 

The function Fi{uj) is strictly positive, while F2{uj) is bounded from below. Further we note 
that the renormalization indicated in (4.14) has been incorporated in the definitions of these 
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functions in such a way each of the individual integrals appearing in (4.17) and (4.18) is 
manifestly finite. 

The functions Fi{uj) and F2{uj) appearing in the action (4.16) are obtained in general 
by numerically performing the integrals in (4.17) and (4.18). However, in the extremal D4- 
brane background (with u-kk = 0) an analytical expression can be given for the renormalized 
worldsheet action (4.16). Using the embedding (2.17) of the D8-D8 pair in the extremal 
background, the renormalized worldsheet action is 



(4.19) 
in the limit that Zoo -^ oo. Thus, the expectation value of (O) takes the form 

L\^/^ / 1 . /7r\ A5 



The exponential dependence on A5/L is precisely that found in [10] coming from the Polyakov 
action. As described above, the overall factor of Nc comes from the Fradkin-Tseytlin contri- 
bution (3.7) to the action. This term also produces the pre-factor of (L/As)^'^. However, one 
should keep in mind that a complete calculation at this order in a' expansion would require 
evaluating the fluctuation determinant on the string worldsheet. Hence one should expect 
this pre-factor to be modified in a complete evaluation at this order. 

Given the general result for the renormalized action (4.16), the expectation value of the 
operator O is given by 

{O) ~ N, (^^y ^ooico)-'/' exp (^^t2^F,{co)+F,ico)^ . (4.21) 

where implicitly we have again taken the limit Zoo -^ 00. The functions Fi{lo) and F2{uj) 
must be determined numerically in the nonextremal background with Ukk 7^ 0. As a check 
of our numerical calculations, the action for the a; = case was determined numerically and 
compared with (4.19). In addition, the Euler number was calculated numerically for each 
a; 7^ embedding and compared with the expected value: x = 1- Iii both cases the numerical 
error, expressed as a fraction of the expected result, was of order of 10^^'' to 10~^^. 

As before, the Polyakov action (i.e., the -Fi term in (4.21)) dominates in the supergravity 
limit and so we focus on this term in the following. The result depends on all three of the 
independent parameters, A5, r^ and L (or rather dimensionless ratios of these parameters) 
— recall that u; is implicitly defined as a function of the ratio — by the relation (4.5). A 
natural approach is to hold the gauge theory parameters constant (by fixing the ratio As/r^) 
and consider the expectation value as a function of L, the separation of the defects. This is 
illustrated in figure 5, where we show log(C') ~ -^ 2w -^i('^) ^^ ^ function of L. The plot 
shows that as L/r4 approaches zero, our result follows the extremal result, tan(7r/16) g:^, 
appearing in the exponential factor in (4.20). This behaviour arises because for L <C r-^, the 
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Figure 5: The contribution to log(C') from tlic rcnormalizcd Polyakov action as a function of L, with 
r4 and A5 fixed. The dashed red curve is the result for the extremal background. The dashed vertical 
line dX L — irr^ corresponds to the defects positioned at antipodal points on the x^ circle. The value 
of log(C') for this configuration is indicated by the dashed green line. The dashed purple line shows 
the value for the closed Wilson line (W) . 

D8/D8 branes do not extend very far into the bulk and so the string stretched between them 
detects no difference between the extremal and non-extremal backgrounds. In terms of the 
gauge theory, this behaviour simply reflects the fact that the chiral symmetry breaking scale 
Uq is much larger than the confining scale C/kk, where supersymmetry breaking takes effect 
in the gauge theory. As L becomes larger, the calculation begins to probe regions of the 
dual spacetime geometry closer to u = Ukk and one sees that the extremal and nonextremal 
behaviours of log(C') begin to deviate around L :^ r^. The expectation value reaches an 
interesting local minimum at L ~ 1.075 r^, where Uq/Ukk — 1-608. Note that the location 
of the minimum is independent of As/r^ and is therefore always visible in the supergravity 
limit. 

When L reaches vrr^, the defects are located at antipodal points on the x'^ circle. This 
corresponds, in the relation (4.5), to the limiting value to = 1. However, we have extended L to 
the region vrr^ < L < 27rr^ in figure 5. Of course, in this regime, the shortest distance between 
the defects on the x^ circle is L = 27rr^ — L but the open Wilson line stretches the longer 
distance L around the circle.^ The embedding profile of the D8-D8 pair is identical to those 
in section 4.1 but with L replacing L. Now in the expectation value, the dual worldsheet 
spans the minimal surface 'outside' of the U-shape formed by the D8-D8 pair. Hence the 
(renormalized) Polyakov action may be calculated as the action of a worldsheet covering the 
entire u-x^ geometry (4.6) minus that for the worldsheet stretched 'inside' of the U-shape. 



*In principle, the following construction could be extended to consider Wilson lines which connect the 
defects after fully winding around the a;* circle some number of times. 
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Figure 6: The contribution to log(C') from the Polyakov action as a function of r^, with L and A5 
fixed. For large r^ the value approaches the result for the extremal background, shown as a dashed 
red line. 

As a result, at this order, we have the relation: {0){L/r^ = 27r — L/r^) x {0){L/r^) = (W) 
where (W) is the expectation value of a closed Wilson line which winds once around the x* 
circle. (Note that we find log(>V) ~ A5/(97rr^).) Figure 5 displays a symmetry about L = irr^ 
which reflects this relation and we may infer that {0){L/r^) approaches zero as L ^ 27rr^. 
Of course, we should add that the five-dimensional gauge theory is defined with a cut-off Uoo 
and so one should not really consider the above results for L,L ^ 1/Uoo- 

An alternative approach to considering the expectation value (O), is to consider it as a 
function of r^ with fixed L and A5. The result is plotted in figure 6 for r^ > L/(j). In the 
decompactification limit, r^ -^ 00, the expectation value again asymptotes to the extremal 
result tan(7r/16) g;;^- Note that, because ui is independent of the ratio A5/L, this plot 
can also be understood as showing the dependence of (O) on the four-dimensional 't Hooft 
coupling, using the relation A4 = A5/(27rr4). 

5. Discussion 

We have examined various aspects of a recent proposal [10, 11, 12] to add quark masses to 
the Sakai-Sugimoto model with nonlocal operators of the form (3.1). The underlying micro- 
scopic field theory is a five-dimensional gauge theory where the chiral quarks are localized on 
separate four-dimensional defects. However, the five-dimensional gauge theory is only defined 
with a cut-off, i.e., new degrees of freedom appear in the far UV. In the dual supergravity 
background, this issue is realized by the running of the dilaton which produces large string 
coupling in the asymptotic region. In section 3, we examined modifications introduced by 
the coupling of the dilaton to the string worldsheet. In particular, we showed that this cou- 
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pling calls for a modification of the renormalization of these operators as in (3.22). The first 
two subtractions, which are linear in the length L, renornialize the Wilson line and are not 
particular to the present open Wilson line calculations. Hence both of these terms, including 
the second one proportional to yju^o, would appear in calculations for closed Wilson lines as 
well. On the other hand, the log subtraction is distinctive of the two end-points of the open 
Wilson line. 

It is interesting to re-express the subtractions in (3.22) in terms of an energy cut-off, 
using the standard dictionary [/qo = 'Uoo/27r£s, 

Iws ' =Iws-LUoo- —^ 77^^ - Q ^OS 7^ • ^-1 

4V2 geff{Uoo) 8 \U,ubJ 

In the third term, QeffiUoo) is the (dimensionless) effective coupling (2.8) of the five-dimensional 
gauge theory evaluated at the cut-off scale C/qo • Hence the a' expansion on the string world- 
sheet produces an expansion in inverse powers of the coupling g^ff from the gauge theory 
perspective, rather than the 1/N^ expansion as produced by a'-corrections to the supergravity 
action — a similar observation was made about the thermal quark diffusion constant in [35]. 
It is interesting that the energy scale Uoo/9eff{Uoo) appearing in the second subtraction is the 
supergravity energy scale associated with u^o [18]. This is a natural energy scale to appear 
here since fluctuations on the worldsheet are contributing at this order [26]. 

In section 4, we explicitly calculated the expectation value of the nonlocal fermion bi- 
linear. This expectation value characterizes the chiral condensate in this holographic model. 
As this holographic construction does not permit the construction of a local fermion bilinear, 
this expectation value is the best order parameter to characterize the spontaneous breaking 
of chiral symmetry. Our explicit calculations yield the result given in (4.21). We note that 
(2.6) and (4.4) can be used to express the pre-factor of the Polyakov term as A5/L ~ OeffiUo), 
up to numerical factors, where QeffiUo) is the effective coupling evaluated at the chiral sym- 
metry breaking scale Uq. The dependence of {O) on L/r^ illustrated in figure 5 describes 
the intricate interplay of the supersymmetry breaking (or confinement) and chiral symmetry 
breaking scales in determining the expectation value. Of course, in the absence of supersym- 
metry breaking, the result in the extremal background (4.20) is independent of L/r^ [10]. As 
figure 5 also illustrates, (O) approaches this supersymmetric result in the limit L/r4 -^ 0. 

Studying the theory at finite temperature in this regime, one finds that the chiral symme- 
try breaking and confinement /deconfinement phase transitions are independent [24, 36]. As 
described in section 2, the chiral symmetry breaking is realized in the gravitational dual by the 
D8-D8 pair joining together to form a smooth U-shaped embedding. The deconfined phase of 
the gauge theory is represented by replacing the supergravity background by a D4 black hole 
[15]. The transition between the low-temperature confining phase and the high-temperature 
deconfined phase occurs when [15] 

-Ideconf ^ Tj • ("-^j 
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Figure 7: For small L, uq is sufficiently large that the D8-D8 pair remains outside of the black hole 
horizon at ut, as shown on the left. For larger L, the branes fall through the horizon, as shown on 
the right. 

In the deconfined phase, if uq is sufficiently large, the tension of the D8-branes can support 
the U-shaped embedding against gravitational attraction of the black hole, which has the 
interpretation that the chiral symmetry remains broken in the deconfined phase [24, 36]. 
Chiral symmetry is not restored until a temperature given by 



1 Lr 



-XSB 



27rr4 L 



0.154 
L 



(5.3) 



where Lc = 0.97r^. Above this temperature, the gravitational attraction becomes sufficiently 
large that the D8-D8 pair are pulled into the horizon (and the embedding is trivial, i.e., 
x^ = constant), as shown in figure 7. The phase structure of the Sakai-Sugimoto model is 
summarized in figure 8. 

In the high temperature phase, where the D8-D8 pair is disconnected, the chiral symme- 
try is restored and so this should be reflected in the expectation value. In particular, beyond 
the phase transition of [24, 36], one should have (O) = 0. In fact, this result does arise 
because with the trivial embedding in the black hole background there is no string world- 
sheet connecting the D8-D8 pair with a single asymptotic boundary. The simplest consistent 
worldsheet would extend through the 'Einstein-Rosen' throat and out to the boundary of 
the second asymptotic region in the black hole geometry. Hence this worldsheet would be 
relevant for a correlator of two operators with the second being in the thermofield double of 
the original gauge theory [37]. 

In a similar way, these expectation values are useful for characterizing the different phases 
in theories with many defects, as discussed, e.g., in [38]. Again, one would find {O) = for 
a Wilson line operator connecting two defects which are not dual to a D8-D8 pair which are 
not joined. 
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Figure 8: Tlic pliase diagram for tlie Sakai-Sugimoto model. The transition from confinement to 
deconfinement occurs at Td — (27rr4)~^. Chiral symmetry is also restored at this temperature if 
L > 0.97 r4. The transitions are independent for L < 0.97 r^, in which case chiral symmetry is 
restored above T^sb — .154/i. 

As observed at the end of section 3.1, {O) oc N^ in accord with the standard large A'^, 
counting. In our calculations, we essentially set the number of flavours to one, however, if 
A^f > 1 one might anticipate the expectation values would also be proportional to N^, again 
reflecting the number of degrees of freedom involved in such a bilinear — e.g., , see [34, 39]. 
However, in the case where Nf > 1, we are implicitly considering the expectation value {O ) 
of an operator with flavour indices / and J for the two fermions. For the smooth embeddings, 
we would have (O^'^) oc 6^^ because consistency requires that the worldsheet start and end on 
the same brane throughout the embedding. A priori, there is no connection between the ip\^ on 
one defect and the V'^ on the other. So our operator reveals this connection as established by 
the chiral symmetry breaking. Tracing over the flavour indices would correspond to implicitly 
summing over the different worldsheets and would produce the factor of A^f mentioned above. 

An alternate approach to understanding chiral symmetry breaking in the Sakai-Sugimoto 
model was considered in [6, 7, 8, 9]. There the key element is the open string tachyon that 
develops between the D8-D8 pair when the (proper) distance separating them is small. Chiral 
symmetry breaking is realized as the condensation of the tachyon, which leads to brane- 
anti-brane annihilation deep in the IR region, producing the smooth embedding in which 
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the D8-D8 pair join. The quark mass and the chiral condensate would be related to the 
asymptotically growing and decaying modes of the tachyon field. This description and the 
approach examined in the present paper both consider the physics of open strings stretched 
between the D8- and D8-branes, so it seems that they must be related. Conceptually, one can 
think of the tachyon analysis as the second-quantized description of the relevant open string 
physics while the worldsheet procedure [10, 11, 12] considered above is the first-quantized 
description of essentially the same physics. Of course, it would be interesting to make this 
connection more precise. This naturally calls for a proper quantization of (open) strings in 
the supergravity background of the D4-brane throat. A more accessible route may be to 
examine the D8-brane embeddings for a nonvanishing quark mass, following the suggestion 
of [10] to include the Polyakov action for the instantonic worldsheet as part of the action for 
the D8-D8 pair . One could then consider the dependence of (O) on niq and compare with 
the results given in [7, 8, 9]. 
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A. Fluctuation Determinant 

In this appendix we will evaluate the contribution from fluctuations for rectangular worldsheet 
discussed in section 3.1. Since a' is a loop-counting parameter on the string worldsheet, the 
one-loop fluctuation determinant will contribute at the same order as the Fradkin-Tseytlin 
term (3.7), considered in the main text. The primary result here is to show that the fluctuation 
determinant contributes no additional UV divergences to the Wilson line calculations, at this 
order. Since the UV behaviour is universal, this result applies for all of the Wilson line 
calculations considered in the present paper. A similar analysis of worldsheet fluctuations 
for a standard closed Wilson line in the D4-brane background has been performed in [40]. 
To begin the calculation, we must expand the worldsheet action to quadratic order about a 
classical solution, X^. Working with the Green-Schwarz formalism,^ this yields 

I = lB + lF + Ighosts, (A.l) 



■''One might question whether or not the Fradkin-Tseytlin term (3.7) is to be added in the worldsheet action 
of the Green-Schwarz string. While the classical action does not couple to the dilaton, this interaction is still 
necessary at the quantum level to preserve the conformal and k symmetry of Green-Schwarz string, just as in 
the bosonic case (and also to have proper effective string coupling dependence of string loops). 
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where 



and 



+ R^^\lsDi^C + hlDiD.^ee]) , (A.2) 

4vr J 

Ighosts = \j <fa^g'^f{g'^^V^e^Vsep - ^R^^^e^ep) . (A.4) 

In the above, i labels a spacetime index which we will split in what follows as a labeling a S'^ 
direction and ^ labeling the remaining transverse directions. The worldsheet directions will 
be labeled by lower Greek indices. The other quantities are specified by [40] 

D^e = d^x\d,e+'cli£!), (A.5) 

A = 5, + \u:,ai>T'''' - ^e'^F^.e^r^^^'^r, ^ a, + M, , (A.6) 

where a, 6, c, d are tangent space indices. The K-symmetry transformation for the GS fermions 
is given by 

5^e = (1 - --=T^pT^^)K . (A.7) 

^v 9 



Note that for simplicity in the following calculations, we will set i? = 1 in the supergravity 
background. We will also drop the dilaton terms in what follows since they are subleading. 
Fermion contributions 

We will use the zehnbeins 

e^" = u^l^dx'' , ^i = l.A (A.S) 

eO = n^/^dt , e^ = u-'^/Uu , (A.9) 

e^ = u^/^dV , e^ = u^'* cos V^X , (A.IO) 

e^ = v}''^ cos ifjsvuxdct)! , e^ = u^' ^ cos ^ cos xc^02 • (A. 11) 

In the following we will sometimes use the notation /Ug = 1, //y = cos V", Ms = cos -0 sin x, A*^ = 
cos V' cos X- The RR field strength in tangent space is given by 

F6789 = - ■ (A. 12) 

u 

In the case at hand, x* = uq and u = ai where (To,(Ti are worldsheet variables. Using the 
pullback metric for the Euclidean worldsheet and noting that \f—g = i\/|ff| we get 

^^^h j "^'"^ (0r"(r«:Po - r,4pi)e - iQ{al'^TJ)i - a~^'^T^,Vo)Q)] . (A.13) 
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Here we have absorbed a' into the fluctuations. Using 

r„Mo - r^Mi = -lu~^^^U + ^^~'/'f r45 , (a.m) 

r„Mi = -V/4f, (A.15) 

o 

r^Mo = lu'/'-lu'/'r, (A.16) 

o o 

with r = r^"^^^ in tangent space we get 

- £!^e(riir4 - ir,)Q + ^^^e(ir + riirr45)e j . (a.i7) 

We will fix K-symmetry with the following: First split Q = 9i + 02- Then choose 

Ts^i = ir4e^ , r^e^ = -iv^e^ , (a.is) 

which leads to after redefining 6^'"^ = u^'^9^''^ and defining d± = ^ia~^''^dQ + di 

ip = ^ j Sa (¥T4d+e^ - Pv^dJ^ + —{¥te^ + Pte^)] . (A.19) 

2tt J \ 4ai ) 

Now choosing the gamma matrices F^ such that 

Fo = ra 1 , Fi = n 1, Fa = rs 7A (A.20) 

with 7A being Euclidean Dirac matrices in 8 dimensions and splitting the ^^'s into two 
Euclidean Majorana-Weyl fermions of opposite chiralities 5", S we get 

If = ^ f d^<J (sd+S -Sd-S + -^{S^S - S^SU . (A.21) 

The squared equations of motion following from the above are: 

[dl + al'^d^al'^di - ^{ido + cj^'^d^) - ^ai)S = , (A.22) 

{dl + al'^dial'^d^ - ^{-ido + a^/^di) - ^ai)S = . (A.23) 

Combining S, S into a worldsheet spinor ^ = {S, S), the equations of motion for the fermions 
then takes on the form 

(di + af^arO, - ^ + ^r^7)^ = . (A.24) 
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Here t^ represents the Pauli matrix ((0, 1), (1, 0)). In first order perturbation theory, the term 
proportional to ri will not contribute. Since S and S combine to form a worldsheet spinor, 
we have 8 massive fermions satisfying the above equation^. 
Boson contributions 

In the second line of (A. 2), we have included the contributions of the Fradkin-Tseytlin 
term (3.7). However, these two terms come with explicit factors of the string length ig which 
reflects the fact that they would only contribute in a two-loop calculation of the fluctuation 
determinant. Therefore we ignore these last two contributions in the following calculation. 
The quadratic order action for the bosons is given by 



A-K 



'' -''"' '-^g + M^(el + « + ^^"^' 



^l + -?C' + ^^l + ^'(^a + -?e + ^d) ) • (A.25) 



Here we have defined ^^4 = u^''^S,x'iAu = u ^''^^u,^fj. = w^'*^C^,la = u^'^^a- Thus the mass 
terms are at 0{l/y/ai). Hence we now have ■^m,^^'*)?^' satisfying 



while ^a satisfy 



(ao^ + afa^.f 5,-^)1 = 0, (A.26) 

ib 



{di + af%ar-%)i = 0, (A.27) 



Ghost contributions 

The ghost action works out to be 



I^kosts = ^/rfV-i^(.1+a?e? + 2^e?+a?(li+a?et+^ei)-24a^/2(e2e^2-eii2)), (A.28) 



with ei = o"i ei and €2 = 0";^ £2. The ghosts satisfy 



32 , 9/2^ -3/2^ I5ai 



id', + al'^d^a-^'^d^ - — i)ei = , (A.29) 

ib 

{di + af/'d,af%-'-^)e, = 0. (A.30) 

Final Result 

The result for the partition function after the above laborious calculation is 

det8/^(a2 + ^3/2g^^3/2g^ _ ^^)detV=^(9f + a^/'d^a~^/'d„ - ^a)det^/\di + ^'3/2^^^9/2^^ _ 15^ 

det6/2(^2 ^ ^3/29^^3/29^ _ 15^)dgt4/2(^2 ^ ^3/25^^3/29^ _ ^^) 

(A.31) 



^Otherwise naively it would appear that there are 8 S"s and 8 S's giving 16 fermions which would lead to 
the wrong counting. 
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The evaluation of the determinant exactly is in general a very hard problem [41]. In this case 
we note that (A. 31) can be rewritten as 

det8/2(a2 + dl- ^)detV2(a2 + dl- §)deti/2(a2 + 9 ^ _ gp ^^ ^^^ 



det6/2(a2 + dl- ^)det4/2(g2 ^ ^2 



b2^ 



where x = --^ and / = 9/4, 51 = 39/4,^2 = 63/4, 61 = 15/4,62 = 3/4. Then each of the 
operators featuring in the determinant can be written as [41] 

so that we can write the determinant as J| A„ where A„ is given by 

4t71 TT 
{-D + — ^^)(/>n = -AnC^-n . (A.34) 

The function 0„ satisfies Dirichlet boundary conditions, namely 0n(a: = 0) = 0, (pnix = Xq) = 
0. Here xq is related to uq through xq = —2/y/uQ. The solution to (A.34) with the Dirichlet 
boundary conditions are known to be Bessel functions ^/xJi{x), ^/xYi{x). Since both it/ 
are allowed, we will choose ^/xJi{x) and ^/xY-l{x) to be the independent solutions. Then 
imposing the boundary condition at x = xq we have 

7/ AiTTi TT 
Jl{u)nXo) = , iX)l_ = {^^f = \n 7^— , (A.35) 



Xo 



L2 



y_,Kxo) = , u^l = (^)2 = A„ - ^^ , (A.36) 

so that a;„'s are related to the zeros ji,n,y-i,n of the Bessel functions. Then the determinant 
can be written using the formula 

00 2 

sinhx = xTT(l + -4-^), (A.37) 

fc=l 



as 



where 



V = VjVy , (A.38) 



sinh^ ^to^ sinhV2 K^^ g.^1^1/2 ^^2^ 

V, = n ^^^L ^:^ci^ ^:^ . (A.39) 

sinh^^kp^sinh^^!^ 
2|xo| 2|xo| 

sinh^ ^^i!^ sinhV2 y_hlf^ sinhV2 ^iff^ 
P, = n '^^ ^ _Ax^L_ . (A.40) 

sinh^ "-'r":^ sinh^ ^~'r":^ 

2|xo| 2|xo| 
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Now we want to get the large n asymptotics of this function. We use the useful identity[42] 
that the large zeros of the Bessel function behave as 

u l^ Au^ -I , ^ ^ 

jy^n cos a - yu,n sm a « (n + - - -jvr - a - — — — ; -. T "^ ' (^-^l) 

2 4 8[(n + I - jJTT - a] 

to get 

n ^ 

-6(n±lT^)-4(n±iT^)), (A.42) 

using which the leading divergence cancels. The subleading terms arise from 0(l/n^) terms 
which lead to a finite result at 0(l/uoo )• The exact formula (A. 38) allows us in princi- 
ple to extract this finite number although we will not attempt it here, as this contribution 
would vanish in the relevant limit Uoo -^ oo. Hence our key result is that in the fluctuation 
determinant (A.31) is in fact precisely 1 in this limit. 
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